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Abstract. 



Q^ ' Asymmetric exclusion processes for particles moving on parallel channels with 

H ■ inhomogeneous coupling are investigated theoretically. Particles interact with hard- 

core exclusion and move in the same direction on both lattices, while transitions 
C^ ' between the channels is allowed at one specific location in the bulk of the system. An 

C/2 . approximate theoretical approach that describes the dynamics in the vertical link and 

horizontal lattice segments exactly but neglects the correlation between the horizontal 

and vertical transport is developed. It allows us to calculate stationary phase diagrams, 

particle currents and densities for symmetric and asymmetric transitions between the 

' O ' channels. It is shown that in the case of the symmetric coupling there are three 

S ■ stationary phases, similarly to the case of single-channel totally asymmetric exclusion 

^ I processes with local inhomogeneity. However, the asymmetric coupling between the 

lattices lead to a very complex phase diagram with ten stationary-state regimes. 

Extensive Monte Carlo computer simulations generally support theoretical predictions, 

K^ ' although simulated stationary-state properties slightly deviate from calculated in the 

^.^ ■ mean-field approximation, suggesting the importance of correlations in the system. 

ly-v , Dynamic properties and phase diagrams are discussed by analyzing constraints on the 

particle currents across the channels. 
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1. Introduction 

In the absence of fundamental framework for non-equilibrium systems, one-dimensional 
models known as asymmetric simple exclusion processes (ASEP), where particles 
interact via an exclusion potential and hop along discrete lattices, have become 
important tools and references for understanding complex phenomena in chemistry, 
physics and biology [H [2]. ASEPs now are playing a role similar to the one the Ising 
model has in the studies of equilibrium systems [3]. Asymmetric exclusion processes 
have been applied successfully to analyze the kinetics of biopolymerization [1], protein 
synthesis [HlEllTj, molecular transport through nanopores and channels [8], the motion 
of motor proteins along cytoskeleton filaments [9l [TOl [11], growth of fungal filaments 
[12] , and different aspects of car traffic processes [131 [H] • 

Many non-equilibrium phenomena can be reasonable well described by a single- 
channel exclusion models. However, the necessity to analyze more realistic complex 
phenomena, that involve the transport along the parallel coupled lattices and existence 
of internal states, e.g. for motor proteins, vehicular traffic and hopping of quantum 
dots, stimulated the development of multi-channel ASEPs [I5l[l6l[l7l[18l[l9l[20l[2l]. In 
these models particle can move along their lattice and they also can transfer between the 
channels, leading to complex dynamic behavior. Theoretical analysis of multi-channel 
exclusion processes indicates that a coupling between channels has a strong effect on 
particle dynamics and stationary-state properties, and it also leads to several unusual 
phenomena, such as localized domain walls [201 [211 [22]. 

Dynamics of many non-equilibrium systems is strongly modified by the presence of 
inhomogeneities and defects. For example, it is known that in the protein synthesis 
the elongation rate of ribosomes is not uniform due to slow codon sites on RNA 
molecules [TJ [23]. The effect of local inhomogeneities have been already investigated 
for different ASEPs on single-chain lattices [6l [2H [251, [26] . However, for multi-channel 
ASEPs mostly homogeneous systems have been considered. The aim of the present 
paper is to investigate two-channel ASEPs with inhomogeneous coupling between the 
lattices. Specifically, we analyze the simplest case when two channels are coupled at 
the single lattice site far away from the boundaries, and two different cases, symmetric 
and asymmetric couplings, are considered. The problem is motivated by the cellular 
transport of motor proteins that move along parallel protofilaments along microtubules 
[27] . Occasionally, motor proteins may jump between the chains, especially in the 
presence of structural defects in the original protofilament. In our study we use a simple 
approximate model for theoretical calculations of stationary-state properties and phase 
diagrams, and we test our predictions by utilizing extensive Monte-Carlo simulations. 

The paper is organized as follows. In section 2 we give a detailed description of the 
model, briefiy present known results for single-chain ASEPs, and develop an approximate 
theory for asymmetric exclusion processes involving two channels connected by a single 
cluster far away from the boundaries where particles can transfer with symmetric or 
asymmetric rates between the lattices. In section 3 we present and discuss Monte-Carlo 
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computer simulation results and compare them with theoretical predictions. Finally, in 
section 4 we summarize and conclude. 

2. Theoretical description of two-channel ASEPs with inhomogeneous 
coupling 

2.1. Model 

We consider identical particles that move on two parallel one-dimensional lattices as 
illustrated in figure 1. Two channels are identical and they have L+l sites. The particles 
interact via hard-core exclusion potential: each site is either empty or occupied by a 
single particle. The dynamics of the system is random-sequential, i.e., at each time step 
we randomly choose a site on one lattice to follow its dynamics. Particles can enter the 
system from the left with the rate < a < 1 if any of the first sites in either lattice 
are empty. Particles that reach the final site in either lattice can leave the system with 
the rate < /? < 1. In the bulk of the system particles move from the site i to the site 
z + 1 on the same lattice with the rate 1 if the next site is available. However, dynamic 
rules are different at the special site L/2 + 1: see figure 1. A particle at this site on 
the first lattice can move vertically to the site L/2 + 1 on the second lattice with the 
rate Wi if the corresponding site is empty, or it can hop to the right empty site with 
the rate (1 — wi). However, if the site L/2 + 1 in the second channel is occupied, the 
particle jumps in the horizontal direction to the right with the rate 1 if the forward site 
is available. The total transition rate out of this defect site is equal to one. Similarly, 
the particle on the site L/2 + 1 on the second lattice can hop vertically with the rate W2 
if the upper site is free, or it can move horizontally with the rate (1 — W2) if the forward 
right site is not occupied. Otherwise, it will move to the right neighboring site with 
the rate 1, assuming again that the forward site is available. When the transition rates 
between the channels are equal {wi = W2) the coupling is symmetric, while for wi 7^ W2 
the coupling between the lattices is asymmetric. Both cases of two-channel ASEPs with 
inhomogeneous coupling are analyzed separately. 

If the inter-channel hopping rates are equal to zero [wi = W2 = 0) then we have 
two independent single-lane ASEPs for which exact solutions for all dynamic properties 
are already known [1]. There are three stationary-state phases. For large entrance and 
exit rates, a > 1/2 and (3 > 1/2, the single-lane ASEP is found to be in a maximal 
current (MC) phase where the stationary current and bulk density are given by 

JmC = T! Pbulk,MC = ^- (1) 

For a < 1/2 and a < f] the dynamics is determined by the entrance process and we 
have a low-density (LD) phase where the particle current and the bulk density are 

Jld = a(l - a), Pbuik,LD = a- (2) 

Finally, for (3 < 1/2 and l3 < a the exiting is a rate-limiting process, and the single-chain 
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Figure 1. Schematic view of two-cliannel ASEP with a coupling at the single site. 
Particles move only along the lattices, except on the link site where they can hop to the 
other lane if the corresponding site is empty. The intcr-channel transitions rates are 
wi and W2- Allowed transitions are shown by arrows. Entrance rates at both lattices 
are equal to a and exit rates are equal to /3. 



ASEP is in a high-density (HD) phase where 

JhD = /5(1 - /5), Pbulk,HD = 1-/9. (3) 

Several approximate theoretical approaches have been developed for investigation 
of multi-channel ASEPs [l5l[T6l[I7l[l8l[T9l[20l[2ll[22]. However, probably the closest 
description of non-equilibrium dynamics, as judged by comparison with Monte Carlo 
computer simulations, is obtained in the mean-field model that have exact treatment 
of the vertical transitions between the channels [Ml [20]. We adopt this strategy to 
analyze the inhomogeneous coupling in two-channel ASEPs. The dynamics at the 
coupling cluster can be specified by introducing functions Pij {i,j = 0,1) that define 
the probability that the vertical cluster is empty (Poo), partially occupied [Pio or Pqi 
if the occupied site is on the lane 1 or 2, respectively], or fully occupied (Pn). These 
probabilities are related via a normalization condition 

Poo + Pio + Poi + Pii = 1. (4) 

In the presented analysis we calculate only bulk densities in all phases, however, our 
analysis can be easily applied to obtain full density profiles utilizing known results [1]. 



2.2. Symmetric coupling 

Let us consider the case of the symmetric coupling when the vertical transitions rates 
are equal. As can be seen from figure 1, the vertical coupling cluster divides the 
two-channel inhomogeneous system into four homogeneous horizontal segments. Each 
of these segments can be viewed as a single-lane ASEP that are connected with the 
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vertical junction cluster, as shown in figure 2. We label left segments in the first and 
second channels as IL and 2L, respectively. The particles enter these segments with 
the rate a and leave them with an effective rate /3e//- Similarly, IR and 2R describe 
the corresponding right segments with entrance and exit rates given by Oe// and /3: see 
figure 2. Because of the symmetry, the probability to have the partially filled vertical 
cluster with the particle on lane 1 or 2 is the same, Pio = Pqi, which modifies the 
normalization condition, 



Pnn + 2Pi 



10 



Pll 



(5) 



Thus dynamics in four segments and in the vertical cluster is treated exactly, but 
correlations between different parts are neglected. 
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Figure 2. Effective scheme for two-channel ASEPs with inhomogeneous symmetric 
couphng. The system consists of four segments coupled by a vertical junction cluster. 
Each segment is viewed as a single-chain ASEP. Allowed transitions are shown by 
arrows. 



It is instructive to analyze first the simplest case of strong symmetric coupling when 
W\ = W2 = ^- It can be shown that it is not possible to have a completely empty vertical 
cluster in the large-time limit, i.e., Pqo = 0. The vertical configuration 00 can only be 
obtained from the partially filled cluster by particles leaving to segments 2L and and 
2R. However, the first particle to reach the junction cluster will hop vertically from 
one defect site to the other constantly until another particle arrives to completely fill 
the cluster. There is no possibility for creating empty vertical cluster. Then from the 
normalization condition ([5]) one obtains, 

1-Pn 



P 



10 



(6) 



Because from the vertical cluster particles can leave to the right part of the system 
only from the completely filled configuration, the effective rate to enter the segments 
IR and 2R can be easily expressed in terms of the properties of the junction cluster, 

«e// = Pll- (7) 
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The effective exit rates from the left segments can also be found using similar arguments: 
the particle can leave the segment IL or 2L to the right if the vertical cluster is in the 
partially filled configuration, leading us to the following expression, 

Peff = PlO = ^^. (8) 

This result has important consequences on the particle dynamics in the system. Since 
Peff < 1/2 the maximal- current phase cannot be realized in the left segments. Symmetry 
and conservation of the particle current also require that 

JiL = J2L = JiR = J2R, (9) 

i.e., fluxes through all segments are equal. Using the notation {x,y), where x and 
y denote the phase of the left and right segments respectively, we conclude that the 
only possible stationary phases for this system are (HD,HD), (LD,HD), (LD,LD) and 
(HD,LD) phases. 

The (HD,HD) phase might exist if the following conditions are satisfied, 

/5<l/2, /3e//<l/2, Peff<a, P<a,jf. (10) 

Because the bulk densities in the left and right segments of this phase are the same, it 
implies that 

Peff = /3- (11) 

Then from Eqs. ([7]) and ([8]) we derive that 

«e// = 1-2/5. (12) 

Substituting this expression into inequalities ( TTOl) we obtain, that (HD,HD) phase is 
realized when 

/?<«, /5<^. (13) 

Similar arguments for (LD,LD) yield that this phase exists for 

a < /3, a < -. (14) 

o 

The conditions for existence of (HD,LD) phase can be written in the following form, 

aeff < 111. Peff < 1/2, aeff < (3, Peff < «• (15) 

In this phase the stationarity of the particle current requires that a^jf = Pejf, and 
utilizing Eqs. ([7]) and (^ it produces 

«e// = Peff = -• (16) 

Combining this result with Eq. (fTS!) we obtain that this phase exists only for 

a > 1/3, p > 1/3. (17) 

It can also be shown that in this mean-field approximation that neglects the correlations 
between the segments the (LD,IID) phase correspond to the phase boundary 

a = p<l/3. (18) 



Inhomogeneous Coupling in Two- Channel Asymmetric Simple Exclusion Processes 7 

In the general case of symmetric coupling between channels, when wi = W2 = w < 
1, the calculations become slightly more involved because the vertical configuration with 
both empty sites now have a nonzero probability, i.e., Pqo 7^ 0. The expressions for the 
effective entrance and exit rates are modified as 

Peff = 1 - Pll - Pio, «e// = Pll + (1 - w)Pio. (19) 

As in the case of strong symmetric coupling, the particle fiuxes through all segments 
are the same, and it is impossible to have a maximal- current phase in the system. If 
MC would exist in one of the segments, than it would have to be found in all segments. 
This would require that aejf > 1/2 and P^ff > 1/2, leading to 

«e// + /5e//>l. (20) 

However, from Eq. ( iT9l) it follows that 

«e// + f^eff = 1 - wPio < 1, (21) 

which contradicts the inequality fl20l) . Thus, our mean- field approach suggests that there 
are four possible stationary phases [(HD,HD), (LD,LD), (HD,LD) and (LD,HD)] in the 
two-channel ASEP with inhomogeneous symmetric coupling. 

Let us define functions pf and P2 as probabilities to find the particle at the last 
sites of left segments IL and 2L respectively. Similarly, pf (z = 1, 2) gives the density at 
the first sites of the right segments 2L and 2R. Because of the symmetry these functions 
are related, 

p^=p^ = p\ pf = p^ = p«. (22) 

The dynamics at the vertical cluster is described at all times by the master equation, 

rIP 

^ = p^(Poi + Pio) - 2Pn(l - P^), (23) 

at 

which at steady state reduces to 

p^Pio = Pii(l-p^). (24) 

The particle current that leaves any of the left segments can be written as 

J = (3,jfp\ (25) 

while the fiux entering any of right segments is given by 

J = a,ff{l-p^). (26) 

Combining Eqs. flTOl) . flMl) . fl25l) and fl26l) . we derive the expression that relates the 
probabilities of filled and partially filled vertical junction clusters, 

(1 - w)Pl + 2PnPio + Pn - Pii = 0. (27) 

Solving this equation and eliminating the unphysical solution leads us to 



P,„ = ^ME^^EMzI^, (28) 

1 — w 
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Then the effective entrance and exit rates can be expressed using Eqs. (1191) as 
«£// = VPii[1 -w{1 - Pii] 



Analyzing the conditions for existence of four possible phase with the help of 
relations (129|) . as was done above in the case of the strong coupling {w = 1), it can 
be shown that (HD,LD) phase exists for 

2-w ^ 2-w , , 

a>- , /3>- . 30 

4 — w 4 — u; 

Bulk density profiles and particle currents in this phase are 

^L _ 2 ^iJ _2-^ . _ . _ 2(2-^) 

(LD,LD) phase is found when 

2 — w 

a < , a < B, (32) 

A-w 

and stationary-state properties are given by 

Pbuik = Pbuik = (^^ Ji = J2 = a{l-a). (33) 

Similarly, (HD,HD) phase exists for 

2 — VI 

p < , p <a. (34) 

A-w 

The bulk density and the particle flux in this phase are equal to 

Pbuik = Pbuik = 1 - /^, Ji = J2 = (3{l- (3). (35) 

The phase (LD,HD) is predicted to be found at the line a = (5 < |5^. However, 
because our theory is a mean-field approach, it can be argued that in the real system it 
does not exist. At the phase transformation line a linear density profile, characteristic 
for first-order phase transitions, is expected [1]. The calculated phase diagram is also 
shown in figure 4. Note that for w = two channels are uncoupled and we recover the 
results for individual single-lattice ASEPs, as expected. In the limit of strong coupling 
{w = 1) the predicted phase diagram is also agrees with the one discussed earlier. 

It is interesting to note that in the case of symmetric coupling the phase diagram for 
two-channel ASEPs is very similar to a phase diagram for the single-chain ASEPs with 
local inhomogeneity [24] . The connection between two systems can be easily understood 
from figure 1. In each lane the hoping rate to the right is the same at all sites except at 
the one that couples the channels. Then the motion along each lattice can be viewed as 
a transport along a single-lane ASEP with local defect that was studied in Ref. |24j . 
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2.3. Asymmetric coupling 

The dynamics in the system of two-channel ASEPs with inhomogeneous couphng 
become very complex for asymmetric vertical transition rates, i.e., when Wi ^ W2-, 
and the symmetry between two lattices no longer exists. Specifically, the densities at 
the last sites of the left segments and densities at the first sites of right segments are 
not equal: pf 7^ pg 5 ^^^ pf 7^ pf . To specify a state of the system we use the notation 
{x,y;w,z) with x, y, w and z corresponding to the phase in the segments IL, IR, 2L 
and 2R respectively, and each of them is able to take the values LD, HD or MC It is 
convenient to view the system with the asymmetric coupling slightly differently from 
the symmetric case, as shown in figure 3. Segments 2L and IR have L/2 sites while IL 
and 2R segments are longer with L/2 + 1 sites. However, because we consider dynamics 
in the thermodynamic limit (L ^ 1), the difference in the sizes does not matter and all 
segments still can be viewed as homogeneous single-chain ASEP with known stationary 
properties. 




lane 1 



Figure 3. Effective scheme for two-clianncl ASEPs with inhomogeneous asymmetric 
coupling. 



To simplify calculations, we consider the case of maximal asymmetry with wi = 1 
and W2 = 0. Our calculations can be easily generalized for general values of vertical 
transition rates. The particle can enter the segment IR only when the vertical cluster is 
fully occupied, and it produces the following expression for the effective entrance rate, 

ae//,i = -Pii- (36) 

However, as can be seen from figure 3, there are two ways to enter the segment 2R: 
from the left segment 2L and from the partially-filled vertical configuration with the 
probability Piq. Then the effective entrance in this case is equal to 



a. 



eff,2 



Pk 



1-pf 



R' 



P^ 



P^^ +Pc 



01- 



(37) 
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Using similar arguments it can be shown that the effective exit rate from the segment 
2L is given by the expression, 

/?e//,l = 1 - ^^, pf = Pii+Pio, (38) 

Pi 

while the effective rate out of 2L segment is 

Peff,2 = 1-pI (39) 

The densities at the last and at the first sites of segments play important role in 
our analysis. They can be easily expressed in terms of the entrance and exit rates from 
known exact solutions for single-chain ASEPs [1]. In the maximal-current phase the 
densities are 

L_ 1 L_ 1 

pf = 1 - Y^. P? = 1 - j-^— (40) 

4tte//,l 4ae//,2 

The corresponding expressions for segments in the low-density regimes are equal to 

]^ a(l — a) j^ a(l — a) 
Pi = —3 'P2 = —3 ' 

Pe//,1 Peff,2 

Pi = Oieff,U P2 = Oieff,2; (41) 

while expressions in the high-density regimes can be written as 

Pi = 1 - /?e//,l, P2 = 1 - /5e//,2, 



1 , Or, i . 



pH^l_ ^^ ^,pK = 1 _ ^iv ^. (42) 

Particle dynamics through vertical junction cluster can be described by a set of 
three independent master equations, 
dPu 



dt 

dPio 
dt 



nnPoi + P2P10 - (1 - Pi)Pu - (1 - n22)Pn; (43) 

(1 - pf - nn)Pn - UuPoi + riu - (1 + nu + P2 )^io; (44) 



dP 

_^ = (l-r222-p^)Pll + (l-p2)A0+P2-(l-^22 + mi + P2)^0i;(45) 

dt 
where un denotes the density at the site before the last one in the segment IL and n22 
is the occupancy of the second site in the segment 2R: see figure 3. At steady-state 
these equations change to the following expressions, 

nuPoi + P2P10 = (1 - Pi)Pii + (1 - n22)Pii; (46) 

(1 - pf - nn)Pn + nn = UuPoi + (1 + ^11 + P2 )Ao; (47) 

(1 - 7222 - P2 )^ll + (1 - P2)^10 + P2 = (1 - ^^22 + Tin + p^)Poi. (48) 
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At large times the stationarity in the system also imphes the following relations 
between the particle currents across different segments: 



tR 



•^1 "T ^2 — ^1 "T "^2 ' "^1 — "^1 



P^ 



10) 



J') — Jo + -P] 



10- 



(49) 



Since each segment can exist in one of three stationary phases (low-density, high-density 
or maximal-current), there are 3^ = 81 possible phases in the system. However, not all 
phases exist due to constraints of steady-state conditions. MC phase cannot exist in the 
segment IR because from Eq. (149|) we have J^ < J^ and 1/4 is the maximal possible 
flux through any segment. Similarly, it is not possible to have MC phase in the segment 
2L since J|' < J^. These arguments already eliminate 27 + 27 — 9 = 45 phases. It is 
also not possible to have at stationary state the phases (HD,LD; x,y) for any x and y. 
For these phase to exist we must have from Eqs. (136|) . (l38l) . (14T|) and (142|) 

Peff,i = 1 - p ^'' and /5e//,i = 1 - Pn - Ao- (50) 

-Til + -T 10 

Comparing these two equations leads to the expression P^q + 2PiiPio = 0, which cannot 
be satisfied at large times. In analogous way it can be shown that phases {x,y; HD,LD) 
are not found at stationary state. Finally, using similar approach we have determined 
that there are only 10 phases that can be supported in the system at t — > cxd. 

The resulting phase diagram is presented in figure 5. Our theoretical calculations 
show that (MC,LD;HD,MC) phase can be found for a > 1/2 and P > 1/2; 
(MC,LD;HD,HD) exists for a > 1/2 and 0.24 < /3 < 1/2; (MC,HD;HD,HD) phase 
is realized for a > 1/2 and 0.22 < f3 < 0.24; (LD,LD;HD,MC) phase exists for 
0.24 < a < 1/2 and (3 > 1/2; and (LD,LD;LD,MC) phase can be found for 
0.19 < a < 0.24 and P > 1/2. The curved boundaries between different phases are 
found by solving the appropriate master equations, and it yields that the curved border 
between (LD,LD;LD,LD) and (LD,LD;LD,HD) phases is given by 



l^ + ^iM 



a)a + -(2a 
4 



l-2a^ -2p + 2P^ + Mf 



-\{M 



l)(3-2/32 + M) 
where the parameter M is defined as 



a(l -a){a-a'^ - (3 + (3'^) 



(51) 



M = VI - 8a + 8a2 + 4/5 - 4/52_ (52) 

The border between the (LD,LD;LD,HD) and (LD,LD;HD,HD) phases is calculated 
from the equation 

1 

2 



2a + M)[(l 



a)a + -{2a 
4 

1 



2a' 



2P + 2p^ + M -1] 



[a 



l){a^ -a + (3-(3^)-- [{M - 1)(1 + 2a^ + 4(3- 4/?^ + M)] . (53) 

The curved boundary between (LD,LD;HD,HD) and (LD,HD;HD,HD) phases is 
computed from 

2a^ + a^ + 2a(3^ - 2a^l3^ + (3^){1 - 2(3 + K) = 



[a 



/3(3 + 2a- 2a^ - 2(3 - K) 



[a — a' 



~(3 + (3^){l + K), 



(54) 
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with 

K = Vl + 4a - 4a2 - gp + S/?^. (55) 

The border between (LD,HD;HD,HD) and (HD,HD;HD,HD) phases is described by 

l + 2a^ + 3(3 + a\l - (3 + (3^) + a{K - 3) = a^ + 3(3^ + K. (56) 

Particle currents and density profiles in each phase can be found by substituting the 
appropriate solutions of master equations into Eqs. ([I]), ([2]) and ([3]). 

3. Monte-Carlo simulations and discussion 

The presented theoretical description of two-channel ASEPs with inhomogeneous 
coupling treats the particle dynamics in horizontal segments and in the vertical junction 
cluster exactly. However, coupling between different parts of the system is viewed in 
a mean-field manner that neglects correlations. In order to test the validity of the 
theoretical analysis we performed extensive Monte Carlo computer simulations. 



1 




T" 



0.8- 



0.6- 



0.4- 



0.2- 



(LD,LD) 



D W=0.1 

V ^=0.5 
O ^=0.9 



(HD,LD) 



99999999 



oooooooooo 



(HD,HD) 



0.4 0.6 

a 



0.8 



Figure 4. Phase diagram for two-cliannel ASEPs witli different symmetric 
inliomogeneous couplings. Symbols correspond to Monte Carlo simulations, while lines 
are theoretical predictions. All phase transformations are first-order transitions with 
density jumps. 



To speed up computer simulations we utilized the so-called BKL algorithm [28] . The 
total number of steps has been varied between 2 x 10^ and 10^. To ensure that our system 
had reached the stationary state we ignored first 3 — 5% of all steps in each simulation. 
Our theoretical results are valid in the thermodynamic limit of L ^ oo, indicating that 
finite-size effects are neglected. In our simulations we used mostly L = 1000 sites lattices, 
but occasionally checked some computations for lattices of sizes L = 100 or L = 10000. 
Computed dynamic properties have been found to be unaffected by lattice size variation, 
confirming absence of finite-size effects. Phase diagram boundaries between different 
phases have been determined by analyzing density profile and particle currents. In the 
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Figure 5. Phase diagram for two-channel ASEPs with asymmetric inhomogeneous 
coupUng {wi = 1 and W2 = 0). Symbols correspond to Monte Carlo simulations, while 
lines are theoretical predictions. Thick lines describe first-order phase transitions, thin 
lines represent continuous phase transformations. 



case where we have transitions between LD and HD phases it is possible to visually 
observe the phase boundary: it happens when the density profile becomes linear. This 
method allowed us to determine the phase boundaries within 0.01 units of the numerical 
values assigned to a and /3. In the case of transitions between an MC phase and LD/HD 
phases, phase borders have been determined by observing the saturation of the particle 
current for specific a and /5 values. This approach, although not as exact as the previous 
one, still allows us to compute the phase transition line to within 0.05 units of the a 
and P values. 

Phase diagram for two-channel ASEPs with inhomogeneous symmetric coupling, 
obtained from Monte Carlo simulations is presented in figure 4, and it can be concluded 
that the approximate theoretical approach agrees well with computer simulations. The 
agreement is excellent for weak coupling between the channels {w <C 1), while increasing 
the probability of vertical transitions {w ~ 1) leads to stronger correlations in the 
system, and theoretical predictions start to deviate from Monte Carlo simulations. More 
information can be obtained from density profiles plotted in figure 6. Our theory gives 
excellent agreement for bulk densities in (HD,HD) and (LD,LD) phases for any coupling 
between the channels: see figures 6a and 6b. This is because the dynamics of the system 
in these cases is governed by entrance or exit processes. However, in (HD,LD) phase 
the theoretical predictions only qualitatively agree with computer simulations, as can 
be seen in figure 6c. This observation can be explained by the fact that in this phase 
the dynamics is governed by processes near the vertical junction cluster, and our theory 
neglects correlations between different segments of the system. 

Stationary-state phase behavior in two-channel ASEPs with inhomogeneous 
asymmetric coupling is much more complex, as shown in figure 5. Monte Carlo 
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Figure 6. Density profiles for two-channel ASEPs with symmetric inhomogeneous 
couplings: a) (LD,LD) phase with w = 1, a ^ 0.25 and /3 = 0.75; b) (HD,HD) 
phase with w = 0.5, a = 0.75 and /? = 0.25; c) (HD,LD) phase with w = 0.9 
and a = l3 = 0.75. Symbols correspond to Monte Carlo simulations, while lines 
are theoretical predictions for bulk densities. 



computer simulations suggest that there are 10 stationary phases, in agreement with 
theoretical predictions. Theoretical calculations for phase boundaries mostly agree with 
computer simulations, although there are deviations for small a and/or small (3. Density 
profiles of six phases for asymmetric coupling are presented in figure 7. The agreement 
between theoretical calculations and computer simulations is excellent in the phases 
where dynamic properties depend on the entrance and exit rates: see figures 7c, 7d and 
7e. However, deviations between the theory and Monte Carlo simulations are observed 
in the phases where dynamics is specified by the processes near the vertical junction 
cluster. Although our theoretical calculations and Monte Carlo computer simulations 
have been performed for the case of strong asymmetry [wi = 1 and W2 = 0), they can 
be easily extended for other asymmetric transition rates, and we expect a similar phase 
diagram and stationary densities and particle currents. 

Our theoretical calculations, supported by Monte Carlo computer simulations, 
indicate that for symmetric coupling between the channels the phase diagram is quite 
simple with only three phases. However, breaking this symmetry leads to a very complex 
dynamics with 10 stationary-state phases. This behavior can be understood by analyzing 
particle fluxes through segments. In the symmetric coupling the currents in all segments 
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Figure 7. Density profiles for two-cliannel ASEPs witli asymmetric inhomogeneous 
couplings (wi = 1 and W2 = 0): a) (MC,LD; HD,MC) phase with a = (3 = 0.75; 
b) (LD,LD; HD,MC) phase with a = 0.4 and /3 = 0.75; c) (LD,LD; LD,MC) phase 
with a = 0.23 and 13 = 0.75; d) (LD,LD; LD,LD) phase with a = 0.1 and /3 = 0.75; 
e) (LD,LD; LD,HD) phase with a = 0.2 and (3 = 0.35; f) (LD,LD; HD,HD) phase 
with a ~ /3 = 0.4. Symbols correspond to Monte Carlo simulations, while lines are 
theoretical predictions for bulk densities. 



are the same. In addition, it was shown that the maximal-current phase cannot exist 
under these conditions. These constraints lead to a simple phase diagram as outlined in 
the figure 4. In contrast, for the asymmetric coupling between the channels the currents 
through all segments are not necessary equal to each other, and MC phase can be found 
under some conditions in some segments. This leads to existence of many steady-state 
phases in the system. 
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The phase diagram for the asymmetric inhomogeneous couphng is almost symmetric 
with respect to a ^^ j3 transformation. This is because the transport of particles from 
the left can be also viewed as a motion of holes in the opposite direction. However, 
as can be clearly seen from figure 5 there is a slight asymmetry in the phase diagram, 
and it also shows up in the dynamic properties. The deviation from the particle-hole 
symmetry can be explained by analyzing the dynamics at the vertical junction for the 
case oi Wi = 1 and W2 = 0, as illustrated in figure 8. The particle at the partially 
filled junction cluster can only move in the vertical direction. At the same time, the 
hole in the partially filled cluster can move vertically and horizontally, if the particle 
from the IL segments enters the cluster. Thus, the corresponding dynamic rules for the 
particles and for the holes in the vertical cluster are not the same, and it produces the 
slight asymmetry in large-time properties of two-channel ASEPs with inhomogeneous 
asymmetric coupling. 
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Figure 8. Dynamic transitions for tlie particle (on the left) and for the hole (on the 
right) at the vertical cluster. Allowed transitions are shown by arrows. 



4. Summary and conclusions 



Two-channel asymmetric simple exclusion processes with inhomogeneous coupling 
between the lattices have been investigated. To better describe the large-time behavior 
of the system a mean-field approach, that treats particle dynamics in four lattice 
segments and in the vertical cluster exactly while neglecting the correlations between 
different parts of the system, is developed. For symmetric transition rates between the 
channels the phase diagram is rather simple and it consists of three phases. This is a 
consequence of constraints on the particle fluxes due to the symmetry. The situation 
is very different for asymmetric coupling between the channels. Here our analysis 
indicates that there are ten stationary phases. It is argued that this complex behavior 
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is realized because of weaker constraints on the particle currents. Our theoretical 
calculations are supported by extensive Monte Carlo computer simulations. Slight 
deviations between theoretical results and computer simulations are only found in the 
phases where correlations around the junction site play important role. Our theoretical 
calculations and computer simulations also point out to asymmetry in the properties 
of two-channel ASEPs with asymmetric coupling with the respect to the particle-hole 
symmetry. It is argued that this asymmetry appears due to different dynamic rules for 
the particles and for the holes at the vertical junction cluster. 

There are several extensions of this work can be proposed. Two-channel ASEPs 
with homogeneous asymmetric coupling show a phase diagram with only seven phases 
and very different stationary properties [22] . It will be interesting to increase the region 
of inhomogeneous coupling from one site to many to see how dynamics of the system 
changes. Another interesting question that will help to understand mechanisms of 
asymmetric exclusion processes is how far coupling clusters should be put in order 
to feel each other. In addition, the analysis can be extended to transport of extended 
particles [7]. 
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